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Abstract 

The production of exotic heavy hadronic particles arises in several models for physics 
beyond the Standard Model. The focus is on R-hadrons, which are stable hadronized gluinos, 
predicted by certain supersymmetric models. Interactions and signatures of single R-hadrons 
are studied with the ATLAS simulation and reconstruction framework. The ATLAS fast 
simulation framework has been extended to include parameterizations for R-hadrons. Based 
on topological and kinematic variables only, the discovery potential of the ATLAS detector 
for R-hadron events produced in pp — > gg is studied for masses below 2 TeV/c^. R-hadrons 
with masses as predicted by standard SUSY scenarios would be discovered already in the 
very early stages of the running of the LHC. The discovery reach of heavy gluinos, predicted 
by for example split supersymmetry models, extends up to at least 1.8 TeV/c^ for three years 
running of the LHC at low luminosity. 



1 Introduction 



Heavy stable hadrons are predicted in several extensions of physics beyond the Standard Model. 
For example, supersymmetry models exist in which the gluino is stable. Models with a stable 
gluino are reviewed in Ref. [1], and include gauge mediated supersymmetry breaking models 
and string motivated super symmetric models. Recently, stable gluinos with masses of the order 
of the TeV scale are predicted in the context of split supersymmetry models [21 [3]. A stable 
gluino would hadronize into heavy (charged and neutral) bound states. These bound states (for 
example gg, gqq, gqqq, qq, qqq) are generically called R-hadrons, where the "R" refers to the 
fact that they can only be stable hadrons if R-parity is conserved [Ij . Several models exist for 
the description of the interactions of stable hadrons in matter [H [5] . The phenomenology of 
stable gluinos has been studied previously in Refs. [H O El [71 [U [9] . 

In addition to supersymmetry, other extensions of the Standard Model have been proposed, 
which predict the existence of new heavy stable hadrons, either due to the presence of a new 
conserved quantum number, or because the decays are suppressed by kinematics or couplings. 
Examples are theories with leptoquarks [TO], theories with universal extra dimensions [llj . the- 
ories with new Standard Model fermions [12] and certain unification models |13j . 

Various searches for stable massive particles have been performed. A comprehensive sum- 
mary may be found in Refs. [141 115j. The negative results of cosmic ray searches and searches in 
matter on earth suggests that these particles would not be strictly stable, so that all cosmological 
bounds can be evaded. In this note, a particle is referred to as stable as long as it does not decay 
in the detector. The mass range from accelerator searches is limited to the center of mass energy. 
In view of the existing mass limits, this study focuses on particles with masses > 100 GeV/c^. 

The detection of heavy non-hadronically interacting charged particles with the ATLAS de- 
tector has earlier been discussed in Refs. [161 117j . In this paper, an overview is presented on 
the interaction and detection of hadronically interacting particles in the ATLAS detector, with 
the focus on R-hadrons. Following the interaction model as proposed in Ref. [5], we discuss 
issues related to triggering, and study R-hadron signatures in the inner detector, calorimeters 
and muon chambers based on full simulation. The discovery potential of the ATLAS detector, 
from an excess of events containing R-hadrons, is investigated for gluinos with masses up to 
2 TeV/c^, using the standard ATLAS detector reconstruction at low luminosity (10 fb~^/year). 
More details about the simulation and analysis of R-hadrons can be found in Ref. |18j . 

The organization of this note is as follows. The production of R-hadrons at the LHC is 
described in Section 2. Section 3 summarizes the event generation for R-hadron and background 
events. Possible trigger signals useful for R-hadron selection are discussed in Section 4. Section 5 
presents a study of the most important signatures from events with R-hadrons in the ATLAS 
detector and describes the final event selection. The resulting mass limits are discussed in 
Section 6. Section 7 addresses specific signatures for R-hadron identification and classification, 
which may also further extend the discovery reach of the ATLAS detector. Finally a summary 
is presented in Section 8. 

2 R-hadron production 

Gluino production is a purely strong process and depends, besides the gluino and squark 
masses, only on the strong coupling constant. A large number of gluinos (roughly 10^ for 
mg=100 GeV/c^, 10^ for mg=l TeV/c^) is expected in one year of LHC running at low luminos- 
ity (10 fb^Vyear). 



1 



2.1 Production channels 



At the LHC, gluinos will be produced via one of the processes pp gg, pp — > qg ^ ggq, or 
pp ^ qq — > ggqq, where the latter two channels are only accessible if the squark is not too heavy. 
For gg production, the leading channels are qq — > gg and gg gg. In this study, we focus on 
the channel gg — > gg, because this process is purely strong and, given the incoming beams and 
energies, depends only the gluino mass and the strong coupling constant. For moderate mg 
(up to ~1400 GeV/c^) this channel is dominating when considering pp gg. The ratio of the 
contribution to gg production from qq to that of gg at the LHC is shown in Fig. [1^ as function 
of rrig and rriq. 
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Figure 1: Aspects of R-hadron production, (a) The ratio of the contribution from qq to that of gg at the 
LHC as function of mg and niq. (b) The distribution of the velocity of the R-hadrons produced at the LHC 
for mg — 300, 900, 1500 GeV/c . Average values are 0.74, 0.64 and 0.56 (c) The distribution of the transverse 
momentum of the R-hadrons produced at the LHC for nig = 300, 900, 1500 GeV/c^ , with average values 188, 427 
and 699 GeV/c. (d) The vector sum of the transverse momenta of the two R-hadrons, with average values of 74, 
129 and 172 GeV/c. For the distributions in (b), (c) and (d) only R-hadrons with < 2.5 were selected. 



2.2 Gluino fragmentation 

Long-lived gluinos hadronize into color singlet bound states. In the hadronization process, 
charged and neutral R-hadrons are produced (gqq, gqqq, gg). For the detection of R-hadrons, 
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it is crucial to know which hadron state is the Ughtest, since that state would dominantly be 
present in the detector (if it is indeed significantly lighter than the other states). The detection 
of neutral hadrons differs significantly from that of charged hadrons. The lightest hadron states 
turn out to be the R-mesons [Sj 119 1 I^O l 121]. which are shown to be nearly mass degenerate, and 
slightly lighter than gluino-gluon states. R-baryons are also expected to be degenerate [SI [20], 
and are roughly 0.3 GeV/c^ heavier than the R-mesons. Thus, there is no preference for an R- 
hadron to be charged or neutral. The probability to form a gg bound state is a free parameter, 
taken to be 0.1 by default [22]. No gg -states are considered in the present study and thus 
R-hadrons are always produced in pairs. For the remaining fraction, the probability to form a 
neutral and charged R-hadron is about half-half. Only a small amount of baryons (~2%) are 
formed due to baryon suppression. In the present study, gluino hadronization is performed using 
dedicated PYTHIA routines [22]. 

2.3 Event topology 

The combination of phase space, matrix elements and parton density functions results in heavy 
particles being produced at a p — p collider like the LHC with high pT values, typically of the 
order of their own mass, i.e. they may be relativistic but their mass is still far from negligible. 
The R-hadrons produced in the channel pp gg events are characterized by two R-hadrons 
produced approximately in a back-to-back configuration in the transverse plane. The likely 
presence of a missing energy signal depends among others on the R-hadron mass, because heavier 
R-hadrons are generally produced with higher momenta than lighter R-hadrons, and thus may 
result possibly in a larger imbalance. The distribution of the velocity and transverse momentum 
of the R-hadrons produced at the LHC via gg gg is shown in Fig. [Tja and c for different 
masses of the gluino. Figure [Tji shows the total transverse momentum of the two R-hadrons 
produced via gg — > gg. 

Due to their slow movement, yet high momentum, and their composite colored structure, with 
one constituent being ultra heavy, R-hadrons posses very distinct features when they traverse 
ordinary matter. The typical signatures for the detection of single heavy stable R-hadrons are 

• High transverse momentum for charged R-hadrons. 

• Large ionization in the tracking system, in case the R-hadrons are charged and slow. 

• A characteristic pattern of energy deposition in the calorimeters. 

• A large time-of- flight, measurable with the muon chambers. 

Isolation may have an importance for the trigger efficiency for R-hadrons, as ATLAS plans 
to employ some form of isolation criteria to control the final trigger rate. Since we expect muons 
to be the dominant background source for R-hadrons, the jet structure around the gluino is 
investigated at event generator level and compared to that of isolated muons and non-isolated 
muons produced via b-quark decay. For this study, 50000 pp ^ Z/7 — > fi^fi~ events and 50000 
pp ^ Z/7 ^ bb events, with b decaying weakly into muons, have been generated, and compared 
to R-hadrons of mass 300 GeV/c^, produced via pp — > gg and pp — > qq ^ ggqq- For the latter 
process the mass of the squark must be larger than the mass of the gluino (LSP) and should 
be below ~1 TeV in standard supersymmetry scenarios. In this example a squark mass of 
800 GeV/c^ is used. In all cases the px of the hard 2^2 process was required to be more than 
50 GeV/c. Figure [2] shows the distribution of the number of particles, the particles density, and 
the total pt of the particles in a cone of size R = y'Arf + Ac/)^ around the muon or R-hadron, 
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Figure 2: A comparison of the jet structure between muons from pp — > Z/7 —> fi'^ jj.^ , muons from pp Z/7 
bb, and 300 GeV/c^ R-hadrons produced via pp gg and pp — > qq — > ggqq. In this example the squark mass is 
set to 800 GeV/c^ Left: mean value of the average number of particles in a cone around the muon direction 
as function of the cone radius R. Middle: mean value of the particle density Ndensity around the muon direction 
as function of the cone radius R. Right: mean value of the transverse momentum, pr, in a cone around the muon 
or R-hadron direction as function of the cone radius R. 

excluding the muon or R-hadron itself (and neutrinos). As can be seen, the typical R-hadron 
will not be perfectly isolated and some loss from an isolation criteria is to be expected. 

3 Event simulation 

3.1 R-hadron event simulation 

There are no models available which propose definite gluino LSP masses following from GUT 
scale assumptions. Thus, the gluino LSP mass at the electroweak scale is a free parameter. 
Below we concentrate entirely on the channel gg — > gg. Since the cross section for this channel 
depends only on the gluino mass, none of the masses of the other supersymmetric particles 
is relevant, as long as they are heavier than the gluino to assure non-decaying gluinos. The 
gluino masses investigated are mg = 100,300,600,900,1100,1300,1500,1700 and 1900 GeV/c^ 
and events were generated with PYTHIA [22j. 

The ATLAS fast simulation framework is used to simulate R-hadrons in the detector and 
perform the physics analysis. The framework is extended to account for R-hadrons, using pa- 
rameterizations of the momentum resolution of inner detector, the response in the calorimeter, 
and the response in the muon spectrometer. The parameterizations come from studying sig- 
natures in the different ATLAS sub-detectors for singly charged R-hadrons of various masses 
(100 GeV/c^ to 1900 GeV/c^), transverse momenta and pseudo-rapidity (r/). These have been 
fully simulated using the ATLAS detector simulation framework based on GEANT3 [23], and 
subsequently reconstructed with the ATLAS C++ reconstruction framework. A detailed dis- 
cussion of the GEANT3 [23] implementation of the model used can be found in Refs. [Sj I18j . 
For the resolution of the inner detector momentum, the parameterization followed quite well 
that of muons in ATLAS [24| , apart from an additional contribution due to multiple scattering 
(proportional to l/(/3p)), especially important for slow R-hadrons, that was included. 

Typically, the momentum resolution, a^l/px), is of the order of 0.005. Details on the 
response of the calorimeters and the muon chambers, forming the base for the parameterizations. 
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are described in Section HI Further information on the parameterizations can be found in 
Ref. pLSj. 

3.2 Background simulation 

QCD processes with their high cross section present a potential background problem for studies 
of hadronic particles. Fortunately, the cross section for QCD events drops rapidly with increasing 
transverse momentum scale in the hard process (pt), and the pT of the hadronized particles is 
further suppressed by the soft fragmentation function of light partons. Moreover, the long 
lifetime of pions and kaons effectively prevents them from decaying into muons before they get 
absorbed in the calorimeter. However, the weak decays of high pT b-quarks and t-quarks will 
produce high pT muons, and thus present a serious background. For the QCD background, 
separate samples have been generated for light quarks (denoted by QCD), b-quarks (denoted by 
bb), and t-quarks (denoted by tt). QCD bb and tt events have been generated in logarithmic 
Pt bins. The number of generated events in each pT bin is assured to be larger than that 
corresponding to 1 fb~^. The only exception to this rule is in the region of very small pT, with 
very large cross sections, but out of several millions of generated events, no events pass the 
trigger cuts. Generating a number of events corresponding to 1 fb~^ is beyond the scope of this 
work, and would only be relevant for a QCD trigger rate study. 

Also, events with single Z and W production as well as events where two gauge bosons 
(denoted by WW/WZ/ZZ) form a potential background source. Such events may result in 
muons, and these events are likely to pass the trigger criteria. The pT value for most muons 
produced from Z and W is of the order of Mz/2, and thus on average smaller than those of 
R-hadrons. Thus, we expect much of this background to be effectively removed by a px cut on 
muons. The generation of Z, W, and diboson events is also done in bins. The width of the bins 
is optimized according to the cross section, and the number of generated events in each pT bin 
is assured to be larger than that corresponding to 1 fb~^. 

4 R-hadron triggers 

Prior to selection, events are required to pass a set of trigger requirements. With their special 
signatures, R-hadrons may fall outside the standard trigger criteria in ATLAS [25^ I26j. In the 
following we discuss the most relevant ATLAS trigger criteria for R-hadron events, as summa- 
rized in Table [TJ The muon triggers turn out to be the most promising, although for very high 
masses, the jet triggers have a sizable contribution. 

The ATLAS trigger efficiencies are estimated by making a pseudo-trigger, which includes 
the trigger menus as defined in Tabled! The pseudo-trigger distinguishes between a 'first level 
trigger', where only calorimeter information and the muon triggers are used, and a 'high level 
trigger', using all possible information. 

4.1 Jet/i?™'''* triggers 

For the jet/S™*^ triggers, it is difficult to estimate the efficiency precisely without dedicated 
full event simulation, because it depends on the details of the specific trigger reconstruction 
algorithms. Two aspects must be kept in mind when studying the interactions of R-hadrons in 
the calorimeter. First of all, an R-hadron passing the ATLAS calorimeters will undergo repeated 
charge flipping in subsequent interactions (~12) through the calorimeter. Hence, the results 
presented here are independent of the initial charge of the R-hadron for the momentum range 
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Menu 


Requirement 


Ex70j70 

jm 

2j350 
^2:200 


^miss > 70 GeV and a jet with Et > 70 GeV. 
Jet Et > 400 GeV. 

Two jets, each of which has Et > 350 GeV. 
^miss > 200 GeV 


2/xlO 


One muon, which has pt > Q GeV/c^ 

Two muons, each of which has pT > 10 GeV/c. 



Table 1: Trigger menus and thresholds used in this study, they represent a subset of tlie full ATLAS trigger 
menus as originally defined in Ref. j25; and later revised in Ref. [26]. 

studied. Second, an R-meson will convert to an R-baryon due to repeated nuclear interactions, 
because this is kinematically favorable [5]. Once the R-hadron is baryon, the nuclear cross 
section is 3/2 larger [5] and thus energy losses due to nuclear scattering increase. 

The energy deposit and the possible stopping of R-hadrons in the calorimeters depend on 
the number of interactions an R-hadron undergoes when passing the calorimeters and the actual 
energy loss per interaction. In Fig. [3^, the number of interactions for an R-hadron traversing 
the ATLAS calorimeters is displayed as a function of r/. The typical energy loss per interaction 
in iron for an R-hadron with a mass of 300 GeV/c^ is shown in Fig. [3)3. As can be seen from this 
plot, the energy loss for an R-hadron tends to level off for very high energies. For an R-hadron 
which punches through the calorimeter and passes the support structure to arrive at the first 
muon-station the number of interactions is about 15% higher. 

The total energy deposit measured in the ATLAS calorimeters is displayed as function of 
rj for various R-hadrons masses and transverse momenta in Fig. HI This figure shows that the 
relative energy loss drops with higher kinetic energy. 

The first level ^et/E^''^^ triggers are based on a quick analysis of the regions of interest 
(ROI's) in the calorimeters and the muon system, while at higher trigger levels this information 
is combined with tracking information [26j . Thus, for R-hadrons, it is crucial to understand the 
interplay between the information accessible at the different trigger levels: for example, it could 
happen that one charged and one neutral R-hadron are produced in the primary interaction. 
Depending on the details of the trigger algorithm and whether a charged R-hadron reaches 
the muon chamber, this could result in a large EJp^^'^ . However, no such imbalance is present 
if considering the calorimeters alone, and it is therefore possible that such an event does not 
pass the first level calorimeter based EJp^^'^ trigger. The same situation may happen for the 
calorimeter jet-based triggers: an R-hadron energy deposit rarely will give rise to a very high 
energetic jet in the calorimeters, but including the possible inner detector track information for 
charged R-hadrons in the final event reconstruction will result in a high pT jet. 

Table [2] gives the final trigger efficiencies, after the first and high level triggers, for the 
jet/S™'^'' trigger menus listed in Tabled! As seen from Tabled an R-hadron event is most 
likely to pass the trigger when the energy deposit in the calorimeters results in a combined 
moderate energy jet and E^^'^^ trigger signal. Due to the moderate energy loss of an R-hadron 
in the calorimeters, combined with the high jet energy trigger threshold, the jet triggers only 
start to contribute for very high R-hadron masses. In particular, the efficiency of the two-jet 
trigger is negligible. The EJp^^'^ trigger is also only effective for higher R-hadron masses. This is 

^Introduced in Ref [25], but replaced in Ref [26] due to the higher (doubled) expected initial luminosity at the 
LHC startup. The local threshold for individual muons in the ATLAS muon-system is still pr > Q GeV/ c. This 
has no impact on the R-hadron signal. 
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Fi gurc 3: (a) The average number of interactions and its spread at the exit of the ATLAS calorimeter system 
as function of rj, for a R-hadron punching through the ATLAS calorimeters. These results are for full simulation 
with the ATLAS detector simulation framework based on GEANT3 and reconstruction with the ATLAS C++ 
reconstruction framework, (b) The average energy loss per interaction in iron for an R-hadron with a mass of 
300 GeV/c^ The lower and upper curves correspond to the energy loss enabling only 2^2 and 2 — > 3 scattering, 
respectively, and the middle histogram represent the combination [5]. 
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Figure 4: Energy loss for punch-though R-hadrons. On the left-hand side the absolute energy loss is shown as 
function of rj for different transverse momenta. The shaded area represents the spread in the energy deposit for 
R-hadrons with mass 300 GeV/c^ and pr of 500 GeV/c. On the right-hand side the relative energy deposition is 
plotted for R-hadron masses of 300, 600 and 1000 GeV/c^ with a transverse momentum of 500 GeV/c. Results are 
obtained by full simulation with the ATLAS detector simulation framework based on GEANT3 and reconstruction 
with the ATLAS C++ reconstruction framework. 

expected, since the typical momentum scales of R-hadron events follow the involved R-hadron 
masses, as discussed in Section 12.31 Finally, it is seen that more than 50% of the R-hadron 
events passing a jet/E'^*^'' trigger also pass a muon trigger. 
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Trigger 


Mass (GeV/c^) 






100 


300 


600 


1100 


1500 


Ex70j70 


Inclusive 
Exclusive 


6% 
1% 


23% 

5% 


36% 
10% 


47% 
18% 


52% 
26% 


jm 


Inclusive 


0% 


1% 


5% 


9% 


11% 


Exclusive 


0% 


0% 


1% 


3% 


5% 


2i350 


Inclusive 
Exclusive 


0% 
0% 


0% 
0% 


0% 
0% 


1% 
0% 


1% 
0% 


Ex200 


Inclusive 


0% 


5% 


13% 


23% 


27% 


Exclusive 


0% 


1% 


2% 


7% 


12% 



Table 2: An estimate of the final jet/ .E™'"* trigger efficiencies, after first and fiigh level trigger, for pp gg events 
at low luminosity. These results are based on the pseudo trigger simulation of R-hadron events in the ATLAS 
detector. The exclusive trigger efficiencies are for events passing a specific jet/i5™"^ trigger and not the muon 
triggers. 

4.2 Muon trigger 

The possibility to use the muon trigger for the detection of slow charged particles, not interacting 
hadronically, has already been investigated for the central region of the ATLAS detector |16j . 
The trigger efficiencies were estimated by requiring a coincidence between the two involved 
trigger stations within a time window of 18 ns. For low luminosity, the efficiency of the trigger 
was found to be as high as 90% for slow particles with /3 = 0.2. However, such slow particles 
would arrive at the muon stations 125 ns after the bunch crossing, or 100 ns after an ordinary 
muon would arrive. Thus, when the muon trigger condition is met, and data-taking started, 
the information of an event which came four bunch crossings later is recorded, and the relevant 
information of the slow particle from the inner detector and calorimeters is probably already lost. 
To which extent the relevant information is available, depends on the velocity of the R-hadron 
and the detailed timing of the signal in the different subsystems. 

In the following, it is assumed that a priori about 75% of all R-hadrons shows up as charged in 
the muon system irrespective of their original charge and type (R-meson or R-baryon) when they 
were formed in the primary interaction. This is due to the many repeated nuclear interactions 
upstream to the muon system, after which an R-hadron independently of it being an R-meson 
or R-baryon will end up as either a charged (3 states) or neutral (1 state) R-baryon at the muon 
chambers. This effect is included in the final trigger estimates presented in Table [5j 

The ATLAS muon trigger has 2 different configurations using different muon stations, de- 
pending on the instantaneous luminosity [24J. For low luminosity the trigger requires a coinci- 
dence between the two barrel RPC's (Resistive Plate Chambers) positioned at approximately 
6.7 (1st RFC) and 7.5 m (2nd RFC) from the interaction point, or the two endcap TGC's 
(Thin Gas Chambers) located at approximately 14 (2nd TGC) and 14.5 m (last TGC). At high 
luminosity, a coincidence is required between different muon stations: the two barrel RFC's 
situated at approximately 6.7 (1st RFC) and 10 m (last RFC), or the two endcap TGC's placed 
at approximately 13 (1st TGC) and 14.5 m (last TGC). The efficiency of the muon trigger is 
studied by requiring in addition that the R-hadron is recorded in the correct event. 

The conditions of an R-hadron firing the muon trigger and being assigned to the right event 
are then the following 

• The R-hadron arriving at the muon stations must be charged and must have r] < 2.5. 
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Figure 5: The difference in time-of-flight ATOF at RPC2 (7.5 m) and TGC3 (14.5 m) between an non-interacting 
particle traveling with velocity /3 in a certain 77-region, and a particle traveling with the speed of light, as function 
of 77. The bunch crossing time Ibc is the maximum difference allowed by the trigger. 

• The R-hadron has a coincidence in the two trigger stations within the temporal gate of 
the trigger stations, which is taken to be 18 ns [IGj for both barrel and endcap trigger 
stations. For example, for the RPC's at 6* = 90°, this would correspond to /3 > 0.5 for the 
high luminosity trigger stations and /? > 0.14 for the low luminosity trigger stations. 

• The delay of a heavy particle compared to that of a muon must be less than 25 ns, so 
that the trigger is fired within the right event. This puts an upper limit on the maximum 
time-of-flight (TOF) of a particle at the trigger station^. In Fig. [5l the difference in time- 
of-flight ATOF between a non-interacting particle with /3 < 1 and a particle with /3 = 1 
at RPC2 (7.5 m) and TGC3 (14.5 m) is shown. The latter stations are the last involved 
stations for a low luminosity muon trigger. Looking at the central detector region (r/ = 0), 
we see that the particle must have (3 > 0.5 in order to reach the last trigger station (7.5 m) 
in time. For a high luminosity trigger making use of the last trigger station in the barrel 
region (10 m), this limit would be /3 > 0.6. In Fig. [6l the delay at the last trigger stations 
for low luminosity is shown for fully simulated R-hadrons of mass 300 GeV/c^ of different 
Pt values as function of rj. Due to the increased energy loss per nuclear scattering, a 
light R-hadron is more delayed than a heavy one with the same starting velocity, when 
measured in the muon system, as shown in Fig. [TOl 

• The deflection angle of the R-hadron in the magnetic fleld must be smaller (pj- greater) 
than that of the corresponding muon (6 or 10 GeV/c). Since R-hadrons are produced with 
high momenta, this condition is always satisfied. 

The l/u6 and the 2/xlO trigger menus are studied by using fully GEANT 3 simulated single R- 
hadrons. The R-hadrons are then required to satisfy the above criteria and the resulting trigger 
efficiencies are shown in Table [3] for pp — > gg. The shift towards lower values of /3 for higher 
R-hadron masses has a clear negative impact on the efficiencies, in contrast to the jet/E™^'* 

^The TOF information in GEANT is not available by default and has been introduced in the ATLAS detector 
simulation framework based on GEANT3 for this purpose. 
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Figure 6: The difference in time of flight ATOF as function of 77 at RPC2 and TGC3 between an R-hadron of 
mass 300 GeV/c^ and a particle traveling with /3 = 1. The displayed momenta correspond to /3-values at 77 = 
of 0.5, 0.6, 0.7, 0.8, 0.9 and 0.98. 



triggers. The full simulation trigger studies presented here does not include the effect a possible 
isolation criteron. Depending on the initial trigger rates in ATLAS, some form of an overall 
isolation criteria can be introduced, although it is not present in the standalone muon system. 
As seen in Section 12. 3| the typical R-hadron is not produced isolated and using a simplified 
pseudo-trigger implementation of the criteria in Ref. [26] leads to a loss of about 50% of the 
events fairly independent of the R-hadron mass. In contrast, an increase of trigger threshold for 
the muon triggers has no impact on the R-hadron signal. 

In the calculations of the efficiencies of the muon trigger, all charged R-hadrons reaching the 
muon chambers are assumed to produce hits in all the trigger planes. Due to support structures 
and the presence of the magnet, it is possible that additional nuclear reactions may take place. 
For high pT muons, the probability to undergo a nuclear interaction is very small, and the 
efficiency of a muon to fire the trigger is very close to 100% |24j . For an R-hadron, however, the 
probability for nuclear interactions is considerably larger. Since the support structures depend 
on T], the probability has been evaluated with full simulation in five different r] regions. The 
results, given in Table [H show that there is indeed a non-zero probability that an R-hadron 
interacts in the muon system and possibly converts into another R-hadron with a different 
charge. This aspect introduces an extra loss of efficiency for the muon triggers. The maximum 
loss, assuming that any interaction in the muon system implies a lost particle, is shown for the 
trigger in Table El The numbers in Table [3] reflect a worst case scenario, as additional events 
where a neutral R-hadron at the entry of the muon-system becomes charged prior to arriving at 
the trigger stations, has not been included. 
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Trigger 


Mass 


Luminosity 


Efficiency (%) 


Maximum loss (%) 




(GeV/c^) 




PP gg 


(Charge flip) 




100 


Low 


61 


3 






High 


58 


9 


l/i6 


300 


Low 


67 


3 






High 


63 


9 




600 


Low 


60 


3 






High 


59 


10 




1100 


Low 


55 


4 






High 


52 


12 




1500 


Low 


45 


5 






High 


43 


13 




100 


Low 


28 


11 






High 


24 


26 


2/xlO 


300 


Low 


33 


12 






High 


29 


28 




600 


Low 


34 


13 






High 


27 


30 




1100 


Low 


21 


13 






High 


17 


31 




1500 


Low 


16 


13 






High 


13 


32 



Table 3: An estimate of the trigger efficiency for pp — > gg events. The fourth column shows the results for fully 
simulated R-hadrons in the ATLAS detector and requiring the conditions given in Section [4. 2 1 The corresponding 
trigger efficiencies at low luminosity obtained from a fast pseudo trigger simulation agree within about 5%. The 
fifth column indicates the potential maximum relative loss of events passing the trigger, due to nuclear interactions 
inside the muon system. 



region 


rj 


5-10 m (total) 


6.7-7.5 m (LL) 


6.7-9.9 m (HL) 




r/(0,0.5) 


30% 


7% 


17% 


barrel: 


77(0.5,1.0) 


40% 


10% 


30% 




77(1.0,1.5) 


45% 


7% 


15% 


region 


77 


7-15 m (total) 


14.1-14.5 m (LL) 


12.9-14.5 m (HL) 


endcap 


77(1.5,2.0) 


60% 


2% 


3% 




77(2.0,2.5) 


75% 


2% 


3% 



Table 4: Fraction of R-hadrons undergoing one or more nuclear interactions in the muon system. The proba- 
bilities are displayed for the entire system, and for the region in between the two trigger stations at low (LL) and 
high (HL) luminosity. 

4.3 Effects of time delays on the R-hadron measurement 

Although it is required that the R-hadrons are in time to be recorded within the right event, the 
event reconstruction may face additional problems as it is not a priori certain whether the full 
detector information of all the separate subsystems can still be used when reconstructing the 
complete event. Below we discuss the time delays in the different ATLAS subdetectors, starting 
from the detector systems closest to the interaction point and ending with the muon system. 
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For the pixels and silicon detectors, the delay in the central region is very small (maximum 1 
and 2 ns for pixels and silicon, respectively), since their location is very close to the vertex. 
Second, for these detectors a hit is based on detecting a certain amount of charge [27]. Slow 
R-hadrons are late, but would be heavily ionizing and certainly result in a hit. Third, we only 
make use of these detectors for momentum measurements, requiring only the hit positions, rather 
than the whole signal shape. Thus, the information of these subdetectors is reliable, if indeed 
we require the R-hadron to arrive at the muon chambers with a delay smaller than 25 ns. 

In the case of the TRT, the delay is slightly larger, maximally (for an R-hadron with velocity 
f3 = 0.5) about 4 ns at the outer layers at 1.1 m. However, this delay is only a small fraction of 
the maximum drift time is of 40 ns, and it would correspond to a drift distance of only 200 //m. 
A good reconstruction of the signal is feasible. 

Heavily delayed R-hadrons imply additional complications for the jet/i?™'^'' triggers. In this 
case the R-hadron signal may not be suitable in time for the calorimeter (trigger) electronics 
to reconstruct the full pulse-shape of the signal, resulting in a lower reconstructed energy. For- 
tunately this is expected to have a small effect for two reasons. Firstly, the delay in arrival at 
the calorimeters is typically small (distance less than 4m), such that it is possible to adjust the 
delay timing to account for this if necessary. Secondly, for first level triggers based on Erp^^^ , 
an R-hadron energy deposit which does not result in a ROI, will contribute to an even higher 
missing energy signal. For full reconstruction of the calorimeters, the time-interval over which 
data taking stretches in case of a positive trigger signal is roughly 125 and 175 ns for LAr and 
Tile calorimeters, respectively. Thus the signal can probably be fully reconstructed, although 
the resolution may be worse due to the time-shift w.r.t. a standard signal, for which the detector 
readout and timing have been optimized. 

For the muon system, the measurement by the individual muon stations will not be affected 
by delays, as the ATLAS muon system will run in a continuous data taking model and with a 
long read out time of about 500 ns [28] . Hence measurement by the individual muon detectors 
are not expected to be affected by time delays below 25 ns. However, the large overall dimension 
of the overall muon system could possibly introduce time incoherence effects when assigning a 
slow particle measurement to the bunch crossing that was triggered on. In a worse case scenario 
of a large delay, the most outer muon chambers could assign a slow traversing charged particle to 
another bunch crossing than the most inner chambers. Although possibly important, studying 
this effect requires a proper time simulation analysis and a realistic description of the muon 
trigger at the time of data-taking, and at the time of writing this is not available in the ATLAS 
simulation. For R-hadrons which are successfully fully measured by the muon detectors and 
assigned to the correct event, the subsequent reconstruction is possibly affected. As seen in 
Section 15. 2^ the time delay will imply a larger of the track fit. It is not expected to be a 
problem to relax the software criteria to include the possibilty to reconstruct particles with /? 
less than one. 

4.4 Trigger rates 

The pseudo trigger used in this study assumes the LHC to be running in the low luminosity 
mode, because the cross sections for gluino production are anyway high, and the low luminosity 
muon trigger has a larger efficiency than the high luminosity trigger. 

In Table [5] the numbers of events expected in 1 fb^^ after triggering are summarized for 
different R-hadron masses and background processes. 
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Sample 


Ntrig(l fb-^) 


Mr =100 


3.6 X lO*" 


Mr =300 


2.0 X 10^ 


Mr =600 


3.6 X 10^ 


Mr =900 


2.3 X 10^ 


Mr =1100 


4.9 X 10^ 


Mr =1300 


1.1 X 10^ 


Mij =1500 


3.2 


Mr =1700 


0.9 


Afij =1900 


0.3 


QCD 


6.5x10*^ 


bb 


4.9 X 10^ 


W 


1.1 X 10^ 


Z 


5.3 X 10^ 


diboson 


1.6 X 10^ 



Table 5: The number of triggered events at low luminosity running for R-hadrons of various 
masses and for the background sources for an integrated luminosity of 1 fb^^. 

5 R-hadron event selection 

A schematic overview of the R-hadron overall behavior in the ATLAS detector is shown in 
Fig. [3 

In the following two selection approaches have been investigated. First, a selection is per- 
formed on the basis of event topology alone. This represents what can be done within the present 
analysis framework, i.e. standard event filter, and which does not require a very fine level of sim- 
ulation of the response of the ATLAS detector. Second, another selection procedure is presented 
requiring a more detailed study of the possibilities to measure the time- of- flight or velocity of 
the R-hadron candidate. The standard ATLAS reconstruction is used to reconstruct events 
containing R-hadrons. Hence, only for initially charged R-hadrons, a high pj- charged track will 
be reconstructed in the inner detector. In contrast, all R-hadrons will produce relatively low 
energetic clusters in the calorimeter and R-hadrons arriving as charged at the muon system will 
be reconstructed by the muon system as a high pT muon assuming a /3 of unity. 

5.1 Event topology selection variables 

Three variables related to the global event topology have been investigated. 

The transverse momentum in the muon chambers: The most energetic track reconstructed by 
the muon system is taken as an R-hadron candidate, and an event is required to have at 
least one candidate with very high transverse momentum (ij?). 

The missing transverse energy, EJp^^'^ : A vector sum of the transverse momenta of all particles, 
calculated using full event information, i.e. also taking into account the inner detector and 
the muon chambers. For charged R-hadron candidates detected in the inner detector, the 
Pt of the track is counted in the vector sum, while for R-hadron candidates with no track 
in the inner detector, the transverse energy as measured in the calorimeter is counted. 
However for an R-hadron candidate detected in the muon chambers but not in the inner 
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Figure 7: Possible manifestations of an R-hadron in ATLAS. 



detector, the pT as measured in the muon chambers is counted. If it is detected in both, 
matching is assumed possible, and only the inner detector is counted in that case. 

The total visible transverse energy of the event, E^"^ : A scalar sum over transverse energies of 
all particles, using the same event information content as EJp^'^^. 

Distributions of the three variables for signal and background after high level trigger cuts are 
shown in Fig. [51 Figure shows the EJp^^^ distribution for R-hadron events with different recon- 
structed detector information. Note the significant difference in the average E^^'^^. For events 
where both R-hadrons have no track or muon information, only the calorimeter information is 
available (like the first level trigger), yielding a moderate E^^^'^. In contrast, events where the 
transverse momentum of only one of the two R-hadron candidates is measured by the inner 
detector or muon system, show a large EJp^^^ , due to the inclusion of the one extra unbalanced 
high pt track. The case where both R-hadron candidates are fully measured, either by the inner 
detector or muon system, has again a moderate EJp^^'^ similar to the first case. 

5.2 Time-of-flight selection variable 

The difference in time-of-flight to the muon chambers, ATOF, between an R-hadron candidate 
and a muon traveling at the speed of light, is a more involved variable, requiring a modified 
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Figure 8: The missing transverse energy and the total visible energy, after high level trigger requirements for 
background (top) and signal (R-hadrons with masses of 300 and 900 GeV/c^). The number of events corresponds 
to an integrated luminosity of 1 fb~^. 




ETmiss('3^^)(°"ly calorimeters) E^missC^^V) g . (QeV) 



Figure 9: The missing transverse energy for R-hadrons produced via gg gg for R-hadron mass of 300 GeV/c^ 
in case only calorimeter information is used (left), in case a neutral and a charged R-hadron are detected (center), 
and in case two charged R-hadrons (right). 



15 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Pstart 



S 1 
^ 0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 










M = 300 GeV/c^ 





















0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

P start 



V 1 T 
^ 0.9- 
0.8- 
0.7- 
0.6 
0.5- 
0.4- 
0.3 
0.2- 
0.1 




M = 600 GeV/c" 



0.2 0.3 0.4 0.5 0.6 0.7 



0.8 0.9 1 

Pstait 



Figure 10: The velocity /3end for difFerent R-hadron masses, as measured at the second layer of MDT's (7.11 
m) in the central barrel region {ri — 0.1), as function of the velocity Pstart at the vertex. The error indicate the 
spread on the mean value. 



track reconstruction in the muon chambers. An estimate of the accuracy of the ATLAS muon 
chambers to measure time-of-flight as well as velocity is given in the following. 

Time-of-flight measurements in ATLAS are possible due to the staggered arrangement of 
the MDT's (Monitoring Drift Tubes) in the muon chambers and an excellent timing precision 
(~1 ns) of the MDT's. 

The multiple repeated measurements of a trajectory from the staggered tube arrangement, 
allow a determination of the velocity of R-hadrons, by using the method of Ref. [17] , which is 
based on minimization w.r.t. /? = u/c of a track. The resulting resolution of the R-hadron 
velocity in the production vertex is found to be 

= 0.067 ±0.029. (1) 



It must be noted that R-hadron energy losses in the calorimeter result in a significantly 
smaller velocity at the muon chambers than that at the primary vertex. As a result of the 
fluctuations in energy loss, the error of the production velocity prediction is increased comparing 
to the velocity measured in the muon system itself. 

In Fig. [To] the velocity of R-hadrons at the second (middle) layer of MDT's (r = 7.11 m), 
fiend-, is shown as function of the Pstart, the velocity at the vertex. A distribution of the velocity 
at the muon chambers is shown in Fig. [11] for an R-hadron and a particle without hadronic 
losses. 

The resolution of the velocity measurement by the muon chambers of R-hadrons given above 
is used to parameterize the transverse momentum measurements by the muon chambers. In the 
following, we also use the ~1 ns timing precision of the ATLAS muon system directly by cutting 
on the measured time-of-flight of R-hadrons compared to that of muons. 



5.3 Selection procedure 

The definition of the signal significance, P, follows that of the ATLAS physics performance 
study [24], 
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Figure 11: The distribution of /3end as measured at the first layer of MDT's (6.92 m) when Patart = 0.8, for a 
strongly interacting particle (solid), and for a particle only interacting electromagnetically (dashed), both of mass 
300 GeV/c^. The spread in f3g„d is much larger for a strongly interacting particle, as expected. 



where S and B are the number of signal and background events, respectively. The usual criterion 
for discovery is a signal significance greater than 5 and at least 10 observed signal events. 
Different cuts apply for different R-hadron mass samples. In Table [U the cuts for all mass 
samples are summarized. The corresponding selected number of events for signal and background 
as well as the S/^/B ratios are listed in Table [71 

These cuts are all related to the global shape and kinematics of the events, i.e. on information 
available at high level trigger. 

Requiring the ATOF with respect to a particle with /3 = 1 to be larger than 3 ns (conser- 
vatively > 3(t) yields a powerful rejection against muons. Thus, for this method the ATOF 
replaces the cuts on missing energy and total transverse energy sum, while the cut on the muon 
transverse momentum is unchanged. For heavy R-hadrons (masses above 1 TeV/c^) the (5 is 
typically 0.8 or smaller and a tighter cut on the time delay could be used, but the limiting factor 
for discovery is the gluino production cross-section. Table [8] shows the corresponding signal 
significances. 



6 Results and discussion 



Using only global event variables, available from standard high level event reconstruction. 
Table [7] shows that R-hadrons can be discovered for masses up to 1400 GeV/c^ for an integrated 
luminosity of 30 fb~^ assuming low luminosity muon trigger configurations. In particular, for low 
masses, a signal significance of 5, the usual criterion for discovery, could be reached after only 
a few days of running. Including more specialized information from a dedicated reconstruction, 
such as time-of- flight, improves the sensitivity and pushes the discovery limit for R-hadron 
masses up to 1700 GeV/c^. The numbers in Tableland [8] represent conservative estimates, as we 
do not include the model-dependent contributions from qq — > gg, which, for high gluino masses, 
can be larger than 100%, see Fig. [1^. For example, in the framework of split supersymmetry 
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Mass 


Cut 


(GeV) 


^ sum 

(GeV) 


Pt 

(GeV/c) 


100 


> 40 


> 200 


> 70 


300 


> 100 


> 250 


> 135 


600 


> 160 


> 450 


> 200 


900 


> 250 


> 600 


> 280 


1100 


> 330 


> 750 


> 300 


1300 


> 380 


> 900 


> 350 


1500 


> 430 


> 1050 


> 430 


1700 


> 470 


> 1200 


> 500 


1900 


> 490 


> 1350 


> 580 



Table 6: Selections for the different R-hadron masses. 

with very high squark masses, the discovery potential would extend up to masses of 1600 GeV/c^ 
if only using global event shape variables, and up to 1800 GeV/c^ if time-of-flight information 
would be included. 

Although the high luminosity muon trigger configuration has not be studied, the reduction 
in muon trigger efficiency is expected to be small, as seen in Table [3l A reduction of less than 
50 GeV/c^ in discovery potential is expected. 

There are a number of uncertainties that may influence the precise discovery limit for R- 
hadrons and in the following we discuss the most important sources. 

The probability to obtain a gluino-gluon state in the hadronization process, Pgg, by default 
taken to be 0.1, is directly related to the ratio of charged to neutral R-mesons produced in 
the hadronization process. The extreme value of 1.0 implies that only neutral R-hadrons are 
produced. Due to the repeated charge flipping of the R-hadrons as they traverse the calorimeters, 
the R-hadron signature in the muon system is independent of the initial charge, hence, a selection 
based on muon-like signatures is unaffected by this parameter. For our selection using missing 
and visible energy, and thus relying on inner detector information, the significance is at most 
reduced by a factor 2 for events generated with Pgg = 1, as the momentum imbalance in 
the tracking system, as well as the total amount of visible energy, are smaller on average. A 
significance loss by a factor 2 lowers the discovery reach by about 100 GeV/c^. Although no 
track is present in the inner detector, the selection is still effective, as the R-hadron can still be 
detected in the muon system and thereby contribute to a total energy signal. Another issue is 
the handling of events with neutral R-hadrons in the inner detector and charged R-hadrons in 
the muon system by the event reconstruction. This uncertainty has not been quantified, because 
a detailed knowledge of the final reconstruction software is needed. But it is not expected to be 
a limiting factor for discovering R-hadrons, as the software can be readapted if necessary. 

As shown in Table HI R-hadrons have a non-zero probability to undergo a nuclear interaction 
in the muon system, and thereby possibly changing their charge. This probability can be as 
high as roughly 50% over the full muon system due to the presence of support structures and 
cryostats. As can be seen from Table [H a large fraction of R-hadron events would actually fire 
the muon trigger, and the uncertainty arises thus mainly from the full track reconstruction in 
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Mass 
(GeV/c^) 


Selected background (1 fb ^) 


Selected 
signal 
(Ifb-i) 


(lfb-1) 


(30 fb-^) 


QCD 


bb 


tt 


W 


Z 


WW/ 

wz/zz 


100 


1.3x10* 


1.1x10^ 


1.2x10'* 


1.7x10* 


8.2x10^ 


4.5x10^ 


8.7x10** 


4.1x10* 


22x10* 


300 


6.2x10^ 


66 


7.3x10^ 


2.0x10^ 


37 


51 


8.6x10* 


1.4x10^ 


7.7x10^ 


600 


40 


6.8 


85 


47 


65 


11 


1.7x10^ 


1.2x10^ 


6.6x10^ 


900 


3.1 


0.6 


11 


13 


2.0 


3.2 


97 


17.1 


93.7 


1100 


1.2 


0.2 


4.0 


5.3 


1.2 


1.3 


19 


5.3 


29.0 


1300 


0.8 


0.1 


1.9 


2.3 


0.8 


0.8 


3.9 


1.5 


8.21 


1500 


0.2 


0.1 


0.8 


1.2 


0.5 


0.5 


1.0 


0.6 


3.29 


1700 


0.1 


0.0 


0.5 


0.7 


0.4 


0.2 


0.25 


0.18 


0.99 


1900 


0.1 


0.0 


0.5 


0.5 


0.3 


0.2 


0.05 


0.04 


0.22 



Table 7: Expected number of signal and background events for the selection presented in TableO 



Mass 


Selected signal (1 fb ^) 


5/75(1 fb"') 


5'/VB(30 fb-') 


100 


1.4x10'^ 


5.9x10^^ 


32x10^ 


300 


1.4x10" 


1.9x10* 


10x10* 


600 


2.8x10^ 


9.8x10^ 


54x10^ 


900 


168 


94.0 


515 


1100 


33 


42.0 


230 


1300 


7.0 


11.8 


65.0 


1500 


1.87 


4.24 


23.2 


1700 


0.44 


1.22 


6.68 


1900 


0.08 


0.24 


1.31 



Table 8: Ratio S/y/B for the different R-hadron masses when replacing the cut on -E^iss 
E'^um by a time-of-flight cut ATOF > 3 ns. 

the muon chambers of such events. The actual loss due to this effect depends on the available 
charge exchange processes for a given initial R-hadron charge and the impact of the lost en- 
ergy in the interaction on the subsequent track reconstruction and acceptance. The worse-case 
scenario is that any interaction inside the full muon system implies that the R-hadrons is re- 
constructed badly and completely lost. In that case, the fractional losses of events containing 
two R-hadrons, due to one or more interactions inside the full muon system (3rd column of Ta- 
bleSj), are approximately 38%, 34% and 31%, for gluino masses of 100 GeV/c^, 300 GeV/c^ and 
600 GeV/c^, respectively. For a reduction of 35% on the number of signal events, the discovery 
reach decreases with about 70 GeV/c^. 

Additional uncertainties related to the generation and simulation of events, such as parton 
distribution functions, higher order corrections and background rate uncertainties have not been 
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studied. Furthermore, the cahbration and final understanding of the detector response may 
influence the above results. No attempt has been made to quantify these kind of uncertainties, 
but given the many and very distinct signatures of R-hadrons, they are expected not to have a 
dramatic impact on the abilities for discovering R-hadrons over a broad mass range. 

7 R-hadron identification from additional detector signatures 

In the situation that an excess of events with slow, yet high pT, tracks in the muon system and/or 
high EJp"^^^ , is observed, it is necessary to study additional signatures to confirm the presence of 
R-hadrons. In addition, the identification of R-hadrons can provide valuable information about 
the structure and the parameters of the underlying theory. The production rate of neutral R- 
hadrons and the fraction of charge hipping can determine whether the heavy spectator is a gluino 
or a squark and measure the probability to form gg bound states. In this section we present 
possibilities to identify and classify R-hadrons by using more detailed detector information such 
as the TRT, the E/p ratio, and shower profiles in the calorimeters. 

These signatures may also help to improve the event selection, but the very high selection 
efficiency and background rejection of a time-of-fiight based selection, makes it difficult to im- 
prove the selection for events with at least one well-reconstructed R-hadron candidate in the 
muon detector. Of the remaining events without a muon trigger, up to about 57% (Table [3|) at 
low luminosity for high R-hadron masses, at most 50% will pass the jet /E'^*'''^ triggers (Table[2]), 
prior to any further event-topology or R-hadron specific selection. As seen in Section [6l a factor 
2 (slightly less than for the case of split supersymmetry) increase in the number of signal events 
changes the discovery reach with at most 100 GeV/c^. 

7.1 Signatures in the TRT detector 

For events with at least one charged R-hadron produced in the hard process, the typically slow 
charged R-hadron will deposit large ionization energy signals in the tracking system. Since the 
TRT has been designed for dE/dx measurements, this subdetector has been investigated in more 
detail. A similar possibility exists for the ATLAS pixel detector j27j . 

The ionization energy deposits of a slow charged particle grow approximately like Hits 
in the ATLAS TRT detector are registered using two discriminator levels: a low threshold level, 
corresponding approximately to 200 eV, and a high threshold level, corresponding to a transition 
radiation photon of about 5 KeV. The time of the low threshold crossing is measured with a 
3.125 ns binning. With these measurements, R-hadron identification can be achieved in two 
ways: by counting the number of straw hits exceeding the high threshold (HT) discriminator 
level i24j or by using the Time-over-Threshold method (ToT) p4l [271129] . 

The number of high threshold hits of R-hadrons, compared to those of electrons and muons, 
has been studied using fully simulated electrons, muons and R-hadrons of masses (100 GeV/c^, 
300 GeV/c^, 600 GeV/c^), generated at various transverse momenta and ry-values. 

In Fig. [T2k the number of straws crossed by a particle is displayed as function of r]. The 
number of crossed straws varies with r/, and so does the the average amount of high threshold 
hits. The TRT detector is therefore divided into five r] regions, and the average amount of high 
threshold hits is studied as function of /?7 of the particle. 

As it can be seen from Fig. [T2b . HT hits can play a role for R-hadron identification only 
for (3^ ^ 1) or /3 < 0.7. On the other hand, the number of HT hits can be used to separate 
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Figure 12: (a) The number of crossed straws as function of R-hadron pseudo-rapidity, (b) Average number of 
HT hits as function of /37 of a particle at central pseudo-rapidity. 
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Figure 13: (a) The relation between the dE/dx of a pion, and its average ToT. (b) A 1-dimensional projection 
of the ToT obtained from testbeam pion data (5 GeV/c) |29) for central values (between 0.8 and 1.2 mm) of the 
distance of closest approach, fitted with a Gaussian. The plot is independent of particle type and depends only 
on the P of the particle. 



fast R-hadrons with /3 > 0.7, when they are minimum ionizing, from very high pT muons 
{pt ^ 200 GeV/c) which start to emit transition radiation photons. 

The Time-over-Threshold method (ToT) |24t [271 122] is based on measuring the time in which 
a signal is above the low threshold discriminator level. In general, the more energy a particle 
deposits, the longer the signal will remain above the threshold, and thus, by measuring the 
time elapsed between a signal exceeding and falling below the threshold, it is possible to extract 
information about the dE/dx: the more ionization energy a particle deposits, the longer the ToT. 
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< ToT > (ns) 

Figure 14: The average ToT for a particle hitting 33 straws (the average over all eta values) for pions of 
momentum 400 GeV/c (solid curve) and for R-hadrons with momentum 400 GeV/c and mass 100 GeV/c^ (or 
/3 = 0.95, dashed-dotted), mass 400 GeV/c^ (or (3 = 0.71, dotted), and mass 600 GeV/c^ (or /3 = 0.55, dashed). 
Note that results only depend on the velocity. 

The TRT offers the possibihty to measure the ToTH. The ToT depends on the hit distance to 
the anode in the straw (the closer to the anode wire the track is located the larger the ToT) and 
on the deposited energy. The ToT averaged over all hit distances in a straw, denoted <ToT>, 
is a function of the dE/dx following the Bethe-Bloch description. In Fig. [T3b the <ToT> is 
displayed as function of the theoretical dE/dx. The linear dependence between the dE/dx and 
<ToT> allows the expected <ToT> of an R-hadron to be obtained assuming that the dE/dx 
for an R-hadron follows Bethe-Bloch dependence, that is, {dE/dx)ji oc (di?/dx)(inip) x (1//?^). 
As is shown in Fig. 113b . the <ToT> has a Gaussian distribution. This plot is obtained from 
testbeam pion data (5 GeV/c) [29j, but is independent of particle type and depends only on the 
P of the particle. 

To estimate the separation power of the ToT method between heavy hadrons and for example 
pions, a simulation experiment is done with 100000 measurements, for different /3 values, an 
average ToT is generated and compared to that of ultrarelativistic pions. The results depend 
only on the velocity and are displayed for pions and R-hadrons with momenta of 400 GeV/c 
in Fig. [HI It is seen that for R-hadrons with velocities (3 = 0.55, /3 = 0.71 and (3 = 0.95 a 
separation of 5.2cr, l.Ocr and 1.6a, respectively, can be obtained. Note that for high /3 values the 
R-hadron becomes less ionizing than an ultra-relativistic pion. 

The number of high threshold hits and the <ToT> in the TRT are the only estimators of 
dE/dx available in ATLAS. These estimators are very highly correlated, so that the separation 
power of the two together is only marginally larger than the separation power of the best single 
estimator. 

^This analysis has been done using test beam data with the TRT electronics developed in 1999. Recently, 
the TRT electronics was changed to provide a more stable drift time measurements independent on the signal 
amplitude. This makes ToT measurement less sensitive to the particle dE/dx. 
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7.2 Matching inner detector and calorimeters 

The low energy deposits from an R-hadron passing the calorimeters (Fig. |4]) makes the com- 
parison of the momentum measurement of the track with energy deposits a potential useful 
separation variable against pions and electrons. A simulation has been made to match inner 
detector momenta with calorimeter deposits of fully GEANT 3 simulated single particles. In 
Fig. [T5t the ratio of the measured energy in the calorimeters and the momentum of the track, 
E/p, is displayed for single R-hadrons, muons, electrons and pions in the barrel region. The 
performance of this method improves for high R-hadron masses due to the lower average /3. A 
quantitative study with fully simulated R-hadron events is required to understand the interplay 
from nearby particles in the subsequent reconstruction. Such a simulation is not yet available. 

7.3 Signatures in the calorimeters 

Due to the charge flipping nature of R-hadron nuclear interactions and the model dependent pro- 
duction rate of neutral R-hadrons, the only certain signature of an R-hadron is that it deposits 
a somewhat moderate amount of energy in the ATLAS calorimeters. Therefore, calorimeter 
measurements could be used in events containing only neutral R-hadrons, or to confirm that the 
particle detected in the muon chambers is a hadron, and not a muon. In addition, the calorime- 
ters can also be used to provide time-of-flight information like the muon system. However, given 
the low energy deposits and the delayed arrival of slow R-Hadrons, the time determination will 
have a larger uncertainty. An additional uncertainty is the propagation time on the shower 
development for an R-hadron, which has not been available in the GEANT3 study. 

Although the precise hadronic energy deposit is subject to many uncertainties [5], a few 
general observations can be made concerning the shower profiles. 
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Figure 15: The ratio E/p for R-hadrons, muons, pions and electrons in the barrel region (r; = 0.1). Singly 
charged R-mesons, muons, pions and electrons are generated and reconstructed. The two extreme cases are 
represented by the muons (solid line) and the electrons (dotted line). 
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Figure 16: Lateral energy profile for muons, electrons, pions and R-hadrons of different masses with pr ~ 100 
and 300 GeV/cas obtained with GEANT 3. 



7.3.1 Longitudinal calorimeter profile 

The longitudinal energy deposit profile for fully GEANT 3 simulated R-hadrons has been studied 
and compared to those of electrons and pions in the different calorimeter compartments of the 
ALTAS detector. 

Differences are due to the fact that the typical energy loss per interaction for R-hadrons 
is small relative to the R-hadrons energy. Contrary to energy losses of electrons or pions, the 
R-hadron energy deposits remain fairly constant along the trajectory of the particle, however it 
should be emphasized that this is only approximately true when nuclear interactions are dom- 
inating; for low P (e.g. when particle is significantly slowed down) the usual electromagnetic 
energy loss is dominating. Since R-hadrons typically are produced with high momenta relative 
to their mass and thus high kinetic energy, the constant energy deposit approximation generally 
holds, and the energy deposit is linearly proportional to the thickness of a calorimeter compart- 
ment. In this aspect R-hadron profiles are similar to those of minimum ionizing particles like 
muons, apart from the absolute scale of the energy deposits. For example, an R-hadron of mass 
300 GeV/c^ and momentum 200 GeV/c, the averaged energy deposits in the tile calorimeter 
compartments amount to 2.5, 6.2 and 3.0 GeV/c^, for corresponding thicknesses of 1.4, 3.9 and 
1.8 interaction lengths, respectively. Note, however, that although the average profiles for muons 
and R-hadrons are similar, the event-by-event variation for R-hadrons, due to charge-flipping 
and fluctuations in the energy loss, are considerably larger than for muons. 

7.3.2 Lateral calorimeter profile 

The lateral energy profile is investigated and compared to that of pions, muons and electrons. 
The width of an R-hadron shower is determined as the fractional energy measured inside a 
wedge as function of its opening angle. The energy is calculated as the sum of the calibrated 
cell energies in the ATLAS calorimeters. 

The results are shown in Fig. [16] for particles generated in the central detector region, with 
an r] value in between and 0.1, and an arbitrary cj) direction. Due to the nuclear interaction 
nature of the R-hadron shower, the shower can be quite broad, depending on the R-hadron 
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kinetic energy. It can be seen that for R-hadrons with large kinetic energy, the shower can be 
as wide as that of a pion (for very high energies even broader), while slow R-hadrons exhibit 
showers which are as narrow as muon showers. However, the full shower shape depends highly 
on the amount of other particles produced in the hadronization process. Although the broadness 
of the shower profile could possibly be used in discriminating R-hadrons from other particles, in 
particular from high px muons, a detailed parametrization and understanding of the event-by- 
event fluctuations of the shower width is left for a future study. 

8 Conclusion 

In this paper we have addressed many aspects important for the detection of heavy hadronically 
interacting particles with exotic and complex signatures in the detector: R-hadrons. Using only 
a few general observables and making use of standard detector techniques with the ATLAS 
detector would allow the discovery of R-hadrons with a mass reach of at least 1400 GeV/c^ 
for an integrated luminosity of 30 fb~^ at low luminosity running. R-hadrons with masses as 
predicted by standard SUSY scenarios (M<600 GeV/c^) could already be discovered at very 
early stages of the running of the LHC accelerator. 

A stand-alone precision muon detector with excellent time-of-flight capabilities, such as the 
ATLAS muon spectrometer, improves significantly the discovery reach and allows a more model- 
independent search. Using only the ATLAS muon spectrometer, this paper has shown that 
R-hadrons produced via gg — gg with masses below 1700 GeV/c^, can be discovered at the 
LHC for an integrated luminosity of 30 fb~^ at low luminosity running. Including R-hadrons 
produced via qq gg, the discovery reach of heavy gluinos, as predicted by for example split 
supersymmetry models, extends up to at least 1.8 TeV/c^. 

At high luminosity the acceptance in the velocity, /?, of the R-hadrons is reduced, due 
to changes in the ATLAS muon trigger configuration, but as the analysis can be performed 
essentially on high-p^ muon like signatures, only a small reduction of the discovery potential is 
expected. 

This paper also discusses additional R-hadron detector signatures which potentially can be 
used to improve the discovery reach even further. 

We conclude that the ATLAS detector has excellent chances of finding R-hadrons with 
masses as predicted by standard supersymmetry scenarios as well as those predicted by split 
supersymmetry models. 
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